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ABSTRACT 

The spectral properties of the adjacency matrix provide a trove of information about the structure and function of complex 
networks. In particular, the largest eigenvalue and its associated principal eigenvector are crucial in the understanding of 
nodes centrality and the unfolding of dynamical processes. Here we show that two distinct types of localization of the principal 
eigenvector may occur in heterogeneous networks. For synthetic networks with degree distribution p(q) ~g~ r , localization 
occurs on the largest hub if y > 5/2; for y< 5/2 a new type of localization arises on a mesoscopic subgraph associated with the 
shell with the largest index in the AT-core decomposition. Similar evidence for the existence of distinct localization modes is 
found in the analysis of real-world networks. Our results open a new perspective on dynamical processes on networks and on 
a recently proposed alternative measure of node centrality based on the non-backtracking matrix. 


Introduction 

An issue of paramount significance regarding the analysis of networked systems is the identification of the most important 
(or central) vertices . 1 The centrality of a vertex may stem from the number of different vertices that can be reached from it, 
from the role it plays in the communication between different parts of the network, or from how closely knit its neighborhood 
is. Following these approaches, different centrality measures have been defined and exploited, such as degree centrality, 
betweenness centrality , 2 or the K -core index and associated /("-core decomposition. ' Among those definitions, one of the most 
relevant is based on the intuitive notion that nodes are central when they are connected to other central nodes. This concept is 
mathematically encoded in the eigenvector centrality 4 (EC) of node i, defined as the component f) of the principal eigenvector 
(PEV) f associated with the largest eigenvalue Aj of the adjacency matrix A tJ . EC is the simplest of a family of centralities 
based on the spectral properties of the adjacency matrix including, among others, Katz’s centrality 5 and PageRank / 1 

Apart from providing relevant information about the network structure , 1 the PEV and associated largest eigenvalue play 
a fundamental role in the theoretical understanding of the behavior of dynamical processes, such as synchronization 7 and 
spreading , 8,9 mediated by complex topologies. Considerable effort has thus been devoted in recent years to the study of the 


spectral properties of heterogeneous networks." 1- ' 1 In this framework, Goltsev et al. 14 (see also 15,16 ) have considered the 
localization of the PEV, i.e., whether its normalization weight is concentrated on a small subset of nodes or not. More in detail, 
let us consider an ensemble of networks of size N, with a PEV /, normalized as a standard Euclidean vector, i.e. jf = 1. 
An eigenvector is localized on a subset V of size Ny if a finite fraction of the normalization weight is concentrated on V 
(Li'ev ff ~ 1)) despite the fact that V is not extensive, i.e., Ny is not proportional to N. This includes the case of localization 

on a finite set of nodes (i.e. Ny independent of N, Ny = 1 in the extreme case of localization on a single node), but also the 
case of localization on a mesoscopic subset of nodes for which Ny ~ with J3 < 1. Otherwise, the eigenvector is instead 
delocalized , and a finite fraction of the nodes Ny ~ N contribute to the normalization weight, implying that their components 
are ~ A -1 / 2 . 

In this context, Goltsev et al. 14 study the localization in power-law distributed networks, with a degree distribution scaling 
as P(q ) ~ q~V, for which the leading eigenvalue Aj is essentially given by the maximum between (q 2 ) / (q) and ^ max , where 
is the largest degree in the network. 9,11 For y > 5/2, where Ai ~ ^/q max , Goltsev et al. 14 find that the PEV becomes 
localized around the hub with degree g max . 14 On the other hand, they argue that, for y < 5/2, when Ai ~ (q 2 ) / (q), the PEV is 
delocalized. These observations are relevant in different contexts. Firstly, they point out a weakness of EC as a measure of 
centrality for heterogeneous power-law networks (y > 5/2), because of the exceedingly large role of the largest hub. 19 On 
the other hand, in the so-called quenched mean-field approach 9, 17 to epidemic spreading on networks, the density of infected 
individuals in the steady state can be related to the properties of the PEV. 14 The localization occurring for large y implies that 
the density of infected individuals in the steady state in those processes might not be an extensive quantity, casting doubts on 
the validity of this theoretical approach and on the actual onset of the endemic infected state. 

Flere we show that the localization properties of the adjacency matrix PEV for heterogeneous (power-law distributed) 
networks are described by a picture much more complex than previously believed. In fact, we provide strong numerical evidence 
that the EC in heterogeneous networks never achieves full delocalization. In the case of uncorrelated synthetic networks with 
a power-law degree distribution, we obtain, by means of a finite-size scaling analysis, that for mild levels of heterogeneity 
(with y > 5/2), the EC is strongly localized on the hubs, as previously argued. For high heterogeneity (y < 5/2), however, we 
point out that the EC, as measured by the components of the PEV, is highly correlated with the corresponding node’s degree. 
This strong correlation results in an effective localization on a mesoscopic subgraph, that can be identified as the shell with the 
largest index in the K -core decomposition of the network. 1 The paper of Goltsev et al . 14 is perfectly correct for what concerns 
the case y > 5/2 but, by only considering the possibility of localization on a finite set of nodes, could not detect the mesoscopic 
localization occurring for y < 5/2. In order to overcome the localization effects intrinsic of the EC, a new centrality measure, 
based on the largest eigenvalue of the Hashimoto, or non-backtracking, matrix, has been recently proposed. 16 We observe that 
this new centrality is not completely free from localization effects. Thus, while it almost coincides with the EC for y < 5/2, 
and for y > 5/2 it avoids the extreme localization around the hubs shown by the EC, it is still localized in this case in some 
mesoscopic subset of nodes, whose characterization calls for further research. The extension of our analysis to the case of 
real world networks is hampered by the fact that usually only one network instance is available, which prevents performing a 
finite-size scaling study. Nevertheless, we numerically argue that also for real networks a twofold scenario holds, in which the 
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PEV is either localized on the hubs, or effectively localized on the maximum K -core of the network. 


Results 

Eigenvector localization and the inverse participation ratio 

A full characterization of an undirected network of size N is given by its adjacency matrix 1 A, whose elements take the value 
Aij = 1 if nodes i and j are connected by an edge, and value A ; / = 0 otherwise. The spectral properties of the adjacency matrix 
are defined by the set of eigenvalues A,, and associated eigenvectors f(A,), i = 1,... ,N, defined by 

Af(A) = Af(A). (1) 

Since the adjacency matrix is symmetric all its eigenvalues are real. The largest of those eigenvalues Ai, is associated with the 
principal eigenvector (PEV) which we denote simply by f. 

The concept of the localization of the PEV f translates in determining whether the value of its normalized components is 
evenly distributed among all nodes in the network, or either it attains a large value on some subset, and is much smaller in 
all the rest. While this concept is quite easy to grasp, assessing it in a single network instance is a delicate issue because any 
quantitative definition involves some degree of arbitrariness. The task becomes however straightforward when ensembles of 
networks of different size can be generated. In such a case, the localization of the eigenvector f associated with the eigenvalue 
A can be precisely assessed by computing the inverse participation ratio (IPR), defined as, 14 - 16 

^A = L/f(A). (2) 

i 

In the absence of any knowledge about the localization support, it is possible to determine whether an eigenvector is localized 
(on some subset in the network) by studying its inverse participation ratio, as a function of the system size N and fitting its 
behavior to a power-law decay of the form 

Y A (N)~N- a . (3) 

If the eigenvector is delocalized, i.e. for f] ~ /V l/2 , the exponent a is equal to 1. An exponent a < 1 is evidence that some 
form of localization is taking place. In the case of extreme localization on a single node, or on a set of nodes with size Ny 
independent of the network size N , the corresponding components of the PEV are finite and this implies K| oc ~ <^(1), i.e., 
a = 0 for N —> Finally, if localization takes place over a subextensive set of nodes of size Ny ~ N&, we expect 

Y A (N)~^N~P, (4) 

leading to a decay exponent a = j3. 

Eigenvector localization in synthetic networks 

We study the localization properties of the PEV computed for synthetic power-law distributed networks of growing size, 
generated using the uncorrelated configuration model (UCM), 19 a modification of the standard configuration model 211 - 21 
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Figure 1 . (a) Inverse participation ratio as a function of the network size for the adjacency matrix of synthetic networks with 
different degree exponent y. For large y, the IPR tends to saturate to a constant value for sufficiently large value of N. For 
y < 5/2, on the other hand, the behavior of the IPR can be fitted to Y^ { (N) ~ N~ a , with a < 1/2. The dashed line represents a 
power-law behavior ~ N 1 , corresponding to a delocalized IPR. (b) Inverse participation ratio as a function of the degree 
exponent y for different network sizes N. The plot confirms the presence of transition in the behavior of the IPR, located in the 
vicinity of y = 5/2. 


designed to avoid degree correlations. 22 In order to explore the presence or absence of localization, we analyze the scaling of 
Y\(N) as function of N as discussed above. In Fig. 1(a) we apply this finite-size scaling analysis to synthetic networks with 
different values of y. In this and the following figures, statistical averages are performed over at least 100 different network 
samples. Error bars are usually smaller than the symbol sizes. In the case of large y we observe an IPR tending to a constant for 
large N, confirming the localization on the hubs predicted by. 14 - 15 The situation is however surprisingly different for y < 5/2. 
Thus, while according to Goltsev et a/., 14 we should expect a delocalized PEV and an IPR decreasing as N a with a = 1, we 
observe instead power-law decays with N, with effective exponents a always smaller than 1/2. The change of behavior of 
the IPR can be further confirmed in Figure 1 (b), where we plot the IPR as a function of the degree exponent y, for different 
values of N. While it is clear that for y > 2.7 the IPR tends to a constant asymptotically, slow crossover effects do not allow to 
draw firm conclusions based on numerics about the precise value of y for which the behavior changes. However, since the 
dependence of the largest eigenvalue on N changes for y = 5/2 11 we expect the transition to take place exactly at y = 5/2: 
simulation results are perfectly compatible with this result. 

The behavior at y < 5/2 can be understood mathematically by observing that the largest eigenvalue in this regime, 
A] = (q 2 ) / (q), l] coincides with the largest eigenvalue of the adjacency matrix in the annealed network approximation. The 
annealed network approximation 2 ' 24 consists in replacing the actual, fixed, adjacency matrix by an average performed over 
degree classes, taking the form 


q'jP{qi\q'j) 

aij ~ NP(qj) 


(5) 


where P(q\q') is the conditional probability that a link from a node of degree q points to a node of degree q. 25 For degree 
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Figure 2. (a) Rescaled scatter plot of f\\N(q 2 )} 172 as a function of q, for a synthetic network with y = 2.1 and size N = 10 6 . 
Data fits the expectation for the PEV in the annealed network approximation, Eq. (7), with only small fluctuations, (b, main) 
Scatter plot of the squared PEV components as a function of the /('-core index for the adjacency matrix of a power-law 
synthetic network with y = 2.1 and size N = 10 6 . (b, inset) Scatter plot of the squared PEV components as a function of the 
degree qj in a synthetic network with y = 3.5 and size N = 10 6 . 


uncorrelated networks, with P{q\q') = qP(q)/(q) 26 we obtain an averaged adjacency matrix 


Cli j — 


Mj 

N (q)' 


( 6 ) 


The matrix an is semi-positive definite and therefore all its eigenvalues are non-negative. IS Then considering that Tr(a 2 ) = 
[Tr(a)] 2 = {(q 2 )/(q)) 2 , where Tr(-) is the trace operator, we have that has a unique non-zero eigenvalue A an = (q 2 ) / (q), 
with associated principal eigenvector fff qj. Applying the normalization condition £,• f 2 = 1, we obtain the normalized form 


f? 


qt 

[.N{q 2 )]V 2 ‘ 


Inserting the expression of ff n into Eq. (2) yields 


(7) 


7 Al W~l/^ (3 “ r)/2 , 


( 8 ) 


that is, a decay with an exponent smaller than 1/2, in agreement with the results in Fig. 1(b). Fig. 2(a) confirms that also 
quenched synthetic networks have PEV components proportional in average to the degree. Notice that Eq. (8) is approximately 
true only in quenched networks for y < 5/2, since the condition leading to it, Eq. (7) fails at y > 5/2, see Fig. 2(b,inset). 

A more physical interpretation of the particular distribution of the PEV in power-law networks with y < 5/2, is that the 
PEV becomes effectively localized on the max(imum) K-coie of the network, defined as the set of nodes with the largest core 
index Km in a K -core decomposition. ■ 27 The K -core decomposition is an iterative procedure to classify vertices of a network 
in layers of increasing density of connections. Starting with the full graph, one removes the vertices with degree <7 = 1, i.e. with 
only one connection. This procedure is repeated until only nodes with degree q >2 are left. The removed nodes constitute 
the K = 1-shell and those remaining compose the K = 2-core. At the next step all vertices with degree q = 2 are removed, 
thus leaving the K = 3-core. The procedure is repeated iteratively. The maximum /C'-core (of index Km) is the set of vertices 
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Figure 3. (a) Total weight Wk m of the PEV on the nodes of the max K -core in synthetic networks as a function of size N. (b) 
Inverse participation ratio as a function of the size of the max K -core Nk m - The dashed line represents a power-law behavior 
~ N k x m . We can see the asymptotic behavior Y\ f ~ N K ^, valid for large network sizes. 

such that one more iteration of the procedure removes all of them. The line of argument leading to this interpretation stems 
from combining the results of Ref ., 14 in which it is proposed that, in epidemic spreading in complex networks , 24 infection 
activity is localized on the PEV, with the observations in Ref ., 29 in which the maximum K -core is identified as a subset of nodes 
sustaining epidemic activity for 7 < 5/2. 

We can see this effective localization on the maximum /('-core in different ways. In the first place, in Fig. 2(b,main) we plot 
the squared components ff of the PEV for all vertices against their corresponding /f-core index. From this Figure we conclude 
that all nodes with the largest f t components belong to the max K -core. The size of this max K -core, Nk m , grows sublinearly as 
a function of the network size as Nk m ~ , 27 However, despite this sublinear growth, a finite fraction of the total PEV 

weight is concentrated on this subset. We check this fact in Fig. 3(a): the total weight of the nodes in the max K- core, 

w Km = e /w 

i\iC K\,, 

tends to a constant in the limit of large network size, implying that more than half of the weight of the normalized PEV 
resides over the max K-corc. Finally, the size dependence of the max K-corc translates, from Eq. (4) in an IPR scaling as 
*A, ~Nk m ” 1 ~ 1/V( 3 “1')/ 2 , in agreement with the result obtained from the degree dependence of the PEV components, ~ qi, 
see Eq. ( 8 ). The relation between IPR and max K-corc size is satisfactorily checked in Fig. 3(b), where we observe it to be 
valid for large network sizes. 

For 7 > 5/2, instead, Figure 2(b,inset) confirms the localization of the PEV around the hub , 14,16 displaying a disproportion¬ 
ately large component on the node with the largest degree. Notice that, irrespective of the value of 7 , with high probability the 
hub belongs to the max K-core. What changes in the two cases is that for 7 > 5/2 the hub alone carries a finite fraction of the 
normalization weight (/) ~ <^(1)) while for 7 < 5/2 it carries a vanishing fraction, and all nodes of the max K-corc must be 
considered to have a finite weight Wk m ■ The behavior for 7 > 3 is clearly evident from Fig. 2(b,inset). In the case 5/2 < 7 < 3, 
the accumulation of a finite weight on the hub takes place for sufficiently large N. This effect is observed in Fig. 4, were we 
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Figure 4. Weight of the PEV as a function of the network size in power-law networks with degree exponent y = 2.8. The 
different functions correspond to: total weight of the nodes in max K- core, Wk m ', total weight in the hub, Wh\ total weight in 
the max k'-core, subtracting the hub, Wk m h ■ 

plot the total weight Wk m of the nodes in the max AT-core, Eq. (9), the total weight in the hub, VV//, and the total weight in the 
max K -core, subtracting the hub, Wk m -h- As we can observe from this Figure, the weight at the hub is small for network sizes 
N < 10 6 , but it then starts to increase, to finally take over, for large network sizes N > 10 7 . 

The non-backtracking centrality 

The observations presented here, together with the arguments provided by Martin et a /., 16 hint that the EC is problematic as a 
useful measure of centrality. For large values of 7 , it is affected by an exceedingly strong localization on the hub, arising as 
a purely topological artifact: the hub is central because its neighbors are central, but those in turn are central only because 
of the hub. For small values of 7 , on the other hand, the observed relation ft ~ q-, indicates that the eigenvector centrality 
provides essentially the same information as the degree centrality. As an attempt to correct the flaws of the EC, Martin 
et al. 16 propose a modified centrality measure, the non-backtracking centrality (NBTC), which is computed in terms of 
the non-backtracking matrix. The Flashimoto, or non-backtracking matrix (NBT ), 16,30,31 is defined as follows: an initially 
undirected network is converted into a directed one by transforming each undirected edge into a pair of directed edges, each 
pointing in opposite directions. If the initial undirected network has E edges, the NBT matrix is a 2 Ex 2 E matrix with rows 
and columns corresponding to directed edges i —> j with value B^j i^ m = 8i^ m { 1 — 8; /), 8ij being the Kronecker symbol. The 
components of the principal eigenvector of the NBT matrix, measure the centrality of vertex i disregarding the contribution 
of vertex j. The NBT centrality of vertex j is given by the sum of these contributions for all neighbors of j: f^ BT = 

The elements of the NBT matrix count the number of non-backtracking walks in a graph and hence remove self-feedback in the 
calculation of node centrality, thus eliminating in principle the artificial topological enhancement of the hub’s centrality. 

As Figure 5 shows, however, the NBTC is not free from localization effects: For all values of 7 the NBTC is not delocalized, 
i.e. Y A[ does not decrease as 1 /N when increasing N. This fact can be understood for 7 < 5/2 in view of the previous results. 
The adjacency matrix PEV is localized on the max K- core, which features many mutual interconnections: the centrality of 
a node is only weakly affected by self-feedback, and removing the contribution of backtracking paths has therefore little 
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Figure 5. Inverse participation ratio as a function of the network size N for the NBTC for power-law synthetic networks with 
different degree exponents 7 . The dashed line has slope —1 indicating delocalization. 

i 

effect. This is confirmed by the scatter plot of the NBTC values / ( XI!T as a function of the corresponding components fo of 
the adjacency matrix PEV, computed for the same synthetic networks. Figure 6 (a). For 7 < 5/2 the two quantities are very 
strongly correlated. For 7 > 5/2 instead. Figure 6 (a) shows that the NBT centrality is truly different and uncorrelated from 
the adjacency matrix EC. However, as Fig. 5 shows, the NBT IPR, computed from the components of the NBTC, decreases 
with the system size N more slowly than N 1 . This is indicative that also in this case a localization occurs on a mesoscopic 
subset, whose size grows sublinearly. Figure 6 (b) shows that this localization is not due to a strong correlation between the 
NBT centrality and the degree of nodes, contrary to what happens for the EC for 7 < 5/2. 

Eigenvector localization in real networks 

For real networks, which have fixed size and do not allow for a finite size scaling analysis, localization is necessarily a more 
blurred concept. The value of K A j gauges how localized the PEV is, but it does not permit to unambiguously declare a network 
localized or not. However, also in this case it is possible to detect, as in synthetic networks, the existence of different localization 
modes. We consider here several real complex networks exhibiting large variations in size, heterogeneity and degree correlations 
(see Methods and Supplemental Material, SM, for details). 

The linear relation between /,■ and the degree q-, is not fulfilled in real networks (see Supplementary Figure SF-1), probably 
due to the presence of nontrivial degree correlations (see SM) which are absent in the synthetic networks. The effective 
localization on the ma x-K core is however still present in some cases. In Fig. 7 we plot for these networks the squared PEV 
component jf as a function of the k'-core index. In some cases (HEP, Movies) all nodes in the max /("-core have a comparable 
and large EC (as in synthetic networks for 7 < 5/2), suggesting localization on the max K- core. In other cases (Internet, 
Amazon) one or a few nodes have a disproportionately large value of jf, hinting at a localization around hubs, as in synthetic 
networks for large 7 . 


8/16 






Figure 6. (a) Scatter plot of the NBTC centralities / ( NHT as a function of the corresponding components of the PEV of the 
adjacency matrix /,, in synthetic uncorrelated networks with a power-law degree distribution. Network size N = 10 6 . (b) 
Rescaled scatter plot of the NBTC centralities ff* BT [N(q 2 )] 1 / 2 as a function of g, for a synthetic network with y = 3.5 and size 
N= 10 6 . 


To clarify the phenomenology we report in Table 1 for each of the real-world networks the values of the leading eigenvalue, 
and the factors (q 2 )/(q) and v /<:/ max . The analysis here is complicated by the presence of degree correlations (see SM), which 
invalidate the direct connection 11 between Ai and the largest between ^/g max and (q 2 )/(<?). 14 However, in some cases (Internet, 
Amazon) the leading eigenvalue is much closer to x /q nVdX than to (q 2 ) / (q): This suggests a localization around the hub and 
matches well with Fig. 7. In others the opposite is true: A] is very far from x /q lmx and relatively close to (q 2 )/(q), hinting 
at a localization on the max K- core, again in agreement with Fig. 7. In other cases (P2P, WWW), values are so close that no 
conclusion can be drawn. 

A further confirmation of this picture is provided by the analysis of the NBT centrality. When localization occurs on hubs 
one expects the elimination of backtracking paths to have a strong impact, as self-feedback effects are tamed. In this case we 
expect the ratio between the IPR for the NBTC and the IPR for the adjacency matrix to be small. On the contrary, when the 
localization occurs on the max k'-core, passing from the adjacency to the NBT matrix would not lead to a big change and we 
expect the ratio to be close to 1. Table 1 confirms this expectation: the IPR ratio is small when the leading eigenvalue A] is 
essentially given by y/q max (localization on hubs) while it is close to 1 when Ai is closer to the (q 2 )/(q) factor (localization on 
the max K- core). A visual representation of these results is provided in Figure 8, where we plot the IPR ratio as a function of 
the ratio between Ai and v /V/ max . As we can see, networks in which the PEV is localized in the max K-core are situated in the 
upper right corner of the panel, while the lower left corner shows the networks with localization occurring on the hubs. 

Discussion 

The properties of the principal eigenvector (PEV), and associated largest eigenvalue, of the adjacency matrix defining a network 
have a notable relevance as characterizing several features of its structure and its effects on the behavior of dynamical processes 
running on top of it. Most important among these features is the role of the components of the PEV as a measure of a node’s 
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Figure 7. Scatter plot of squared PEV components of the adjacency matrix of the real-world networks as a function of the 
/C-core index. 


importance, the so-called eigenvector centrality. One of the properties of the PEV that has recently attracted the interest of 
the statistical physics community is its localization. In the case of networks with a power-law degree distribution P(q) ~ q 
initial research on this subject 14,16 suggested that, for y > 5/2, the PEV is localized on the nodes with largest degree. On the 
other hand, for y < 5/2, the PEV should be delocalized. 

In this paper we have shown that eigenvector localization in heterogeneous networks is described by a more complex 
picture. Thus, we present evidence that for all power-law distributed networks the PEV is always localized to some extent. In 
the case of synthetic power-law distributed networks, we observe that, while for mildly heterogeneous networks with y > 5/2 
the PEV is indeed localized on the nodes with maximum degree (the hubs), in the case of high heterogeneity, with y < 5/2, the 
PEV shows a peculiar form of localization, its components /, being proportional to the node’s degree, f) ~ qj. This particular 
proportionality induces an effective localization on the maximum K -core of the network, defined as the core of maximum index 
in a /('-core decomposition. This max K-corc concentrates a finite fraction of the normalized weight of the PEV, despite the 
fact that the size of the max K-corc is sublinear with the network size. In the case of real world networks, the elucidation of 
the PEV localization is not so clearcut. We however provide evidence for an analogous scenario as that observed in synthetic 
networks, where the nature of the localization of the PEV is ruled by its associated largest eigenvalue A|: When Ai is close to 
the mean-field value (q 2 )/(q), localization on the max K -core is expected. On the other hand, when the largest eigenvalue is 
close to y/q max , localization takes place on the hubs. 

The results presented here give a new perspective on complex topologies from several perspectives. Firstly, it is common 
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Figure 8. Ratio between the NBTC IPR and the IPR of the adjacency matrix as a function of the ratio between the largest 
eigenvalue and the square root of the maximum degree, for the real networks considered. The symbol codes are: square for 
localization on the max /f-core; circle for localization on the hub; triangle up for networks in which x /q„y dX is very close to 
(q 2 ) / (q), so no conclusion can be drawn; triangle down for the rest of networks. 

knowledge that networks with 7 > 3 are fundamentally different from those with 7 < 3 (scale-free networks) because the 
divergence of the second moment of the degree distribution has a series of crucial effects. A tacit corollary is that networks 
with 2 < 7 < 3 have essentially the same properties. Our paper, together with other recent results , 14 points out that networks 
with exponent 7 < 5/2 are in many respects qualitatively different from those with 7 > 5/2. Secondly, our results point out 
the weakness of eigenvector centrality as a measure of centrality for power-law networks. Indeed, for 7 < 5/2, eigenvector 
centrality does not provide more information than degree centrality, while for 7 > 5/2 the eigenvector localization on the 
hubs arises as a purely topological artifact. Alternative measures of centrality, based on the Hashimoto non-backtracking 
matrix 16,30,31 are also not free from localization effects. Finally, from a dynamical point of view, largest eigenvalues and the 
associated eigenvectors are crucially related to the properties of processes on networks 7,14,13 and their localization effects 
should be taken properly into account when developing theories relying on the structure of the adjacency matrix. 

The localization properties described here call for a revision of our present understanding of heterogeneous topologies. 
Other networks properties, such as degree correlations, clustering or the presence of a community structure, might play a role 
in the localization of the PEV. The clarification of these effects, as well as the understanding of the nature of the mesoscopic 
subgraph on which the NBTC is localized for 7 > 5/2, are still open questions, calling for further scientific effort. 

Methods 

Real networks analyzed 

We consider in our analysis the following real networks datasets: 

• HEP: Collaboration network between authors of papers submitted to the High Energy Physics section of the online 
preprint server arXiv. Each node is a scientist. Two scientists are connected by an edge if they have coauthored a 
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preprint. 33 

• Slashdot: User network of the Slashdot technology news website. Nodes represent users, which can tag each other as 
friends or foes. An edge represents the presence of a tagging between two users. l4 

• Amazon: Co-purchasing network from the online store Amazon. Nodes represent products, which are joined by edges if 
they are frequently purchased together. 

• Internet: Internet map at the Autonomous System level, collected at the Oregon route server. Vertices represent 
autonomous systems (aggregations of Internet routers under the same administrative policy), while edges represent the 
existence of border gateway protocol (BGP) peer connections between the corresponding autonomous systems. 36 

• Email: Enron email communication network. Nodes represent email addresses. An edge joins two addresses if they have 
exchanged at least one email. 14 

• P2P: Gnutella peer-to-peer file sharing network. Nodes represent hosts in the Gnutella system. An edge stands for a 
connection between two Gnutella hosts. l3 

• Movies: Network of movie actor collaborations obtained from the Internet Movie Database (IMDB). Each vertex 
represents an actor. Two actors are joined by an edge if they have co-starred at least one movie. 17 

• WWW: Notre Dame web graph. Nodes represent web pages from University of Notre Dame. Edges indicate the presence 
of a hyperlink pointing from one page to another. lS 

• PGP: Social network defined by the users of the pretty-good-privacy (PGP) encryption algorithm for secure information 
exchange. Vertices represent users of the PGP algorithm. An edge between two vertices indicates that each user has 
signed the encryption key of the other. 39 

Some of this networks are actually directed. We have symmetrized them, rendering them undirected, to perform our 
analyses. 
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Network 

N 

(?) 

(? 2 )/(?> 

\J Qmax 

Ai 

r(At) 

Y(A N l BT ) 

IPR Ratio 

HEP 

12006 

19.74 

129.94 

22.16 

244.93 

0.003890 

0.003887 

0.9993 

Slashdot 

82168 

12.27 

149.71 

50.52 

134.63 

0.002174 

0.002006 

0.9228 

Amazon 

403394 

12.11 

30.55 

52.46 

57.15 

0.089122 

0.005423 

0.0608 

Internet 

10790 

4.16 

259.46 

48.34 

59.58 

0.066138 

0.015783 

0.2386 

Email 

36692 

10.02 

140.08 

37.19 

118.42 

0.003790 

0.003446 

0.9091 

P2P 

62586 

4.73 

11.60 

9.75 

13.18 

0.000921 

0.000592 

0.6429 

Movies 

81860 

89.53 

594.92 

61.55 

817.36 

0.000640 

0.000638 

0.9966 

WWW 

325729 

6.69 

280.68 

103.54 

184.93 

0.022726 

0.008357 

0.3677 

PGP 

10680 

4.55 

18.88 

14.32 

42.44 

0.016622 

0.015989 

0.9619 


Table 1 . Relevant metrics for the various real-world networks with and the measured value of the IPR ratio between Y(A^ BT ) 
and Y (Ai). Size and other information on the networks are provided in the Supplementary Information. 
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